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Abstract 

Travel time perturbations of adiabatic normal modes due to an internal tide and 
internal mode field in the Barents Sea are examined. A formalism for the travel time 
perturbation due to a change in sound speed is presented. Internal tide and internal 
wave amplitude spectra are calculated from Brancker temperature loggers which were 
deployed on moorings in the Barents Sea during the August 1992 Barents Sea Polar 
Front Experiment. In particular, the first three internal wave mode amplitudes are 
estimated from the four Brancker temperature loggers on the southwest mooring of 
the array. Modal perturbations in acoustic pulse travel time and the travel time 
covariance are calculated and compared for consistency to a simple ray model. These 
perturbations are small for the modal arrivals that the vertical acoustic array which 
was deployed is expected to resolve. The third internal wave mode has the largest 
impact on the acoustic arrivals, per unit amplitude, but the first internal wave mode 
dominates the scattering due to having a much larger amplitude overall. 
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Chapter 1 



Introduction 



Tomography comes from the Greek word tomos meaning slice or section. Ocean 
acoustic tomography is the science of taking acoustic slices or sections of the ocean 
sound speed, temperature, and current fields. In this thesis we will deal with the 
application of tomographic techniques to the study of internal waves and internal 
tides. 



1.1 Motivation 

The ocean covers over two-thirds of the earth’s surface. With its tremendous heat 
capacity and enormous currents, the ocean exerts a significant influence on the world’s 
climate and meteorology. This influence is not well understood, in part due to the 
difficulty in observing the ocean volume, and also in part due to the ocean’s meso- 
scale and inter-annual variability. Remote sensing using electromagnetic energy in 
the forms of light and radar, which has allowed man to peer billions of light-years 
into outer space, encounters in the ocean an opaque medium. The average depth of 
the ocean is 4000 m, but blue-green light and radar can only penetrate to depths 
of 100 m and the surface, respectively. What knowledge we have about the ocean 
depths was obtained mainly through in-situ measurements conducted by a myriad 
of research vessels performing tens of thousands of casts of Nansen bottles, Niskin 
bottles, and CTD’s. 
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Though bottle and CTD casts are important and still continue today, they are 
time and labor intensive. Not only are they costly, but the slow sampling rate makes 
them susceptible to aliasing of many oceanic processes. As new methods have been 
sought out, the use of acoustics to probe the depths of the ocean has become an 
attractive, alternate technology in oceanic research. Acoustics has allowed man to 
determine accurately the topography of the ocean floor, to track marine life, to profile 
the ocean currents, and to perform many other useful tasks. One of the more recent 
research efforts in acoustics has been that of ocean acoustic tomography. The idea 
behind tomography is to take many slices or sections of an object at different angles to 
produce a three dimensional picture of the subject. The theory was mathematically 
proven by Radon in 1917 [cited in Munk and Wunsch , 1979]. Tomography has been 
widely used in medicine (CAT scan stands for Computer Assisted Tomography scan), 
and in seismology. 

The basic principle behind ocean acoustic tomography is that by measuring the 
long-term, temporal variation in acoustic pulse travel times between a source and 
a receiver, due to the sound speed field variation from effects such as fronts and 
eddies, one can image these features acoustically, and thus study them. This temporal 
variation had been noted in several early experiments [see Hamilton , 1977; Steinberg 
and Birdsall, 1966 (both cited in Munk and Wunsch , 1979)], but was often considered 
noise. Worcester [1977] conducted one of the first experiments in ocean acoustic 
tomography, when he measured reciprocal acoustic travel times between two ships 25 
km apart, in order to investigate the feasibility of using the travel time differences 
to monitor mid-ocean currents. Munk and Wunsch [1979] proposed expanding the 
range between source and receiver to 1000 km to monitor mesoscale ocean features. 
Since then, acoustic tomography has been applied to both small scale and large scale 
oceanographic eddies, as well as to shallow water and deep water environments. 
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1.2 Thesis Objectives 



In this work, we will first calculate internal tide and internal wave spectra and am- 
plitudes from temperature measurements taken at moorings at several depths in the 
Barents Sea during the August 1992 Barents Sea Polar Front (BSPF) Experiment. 
Using this oceanographic data, the perturbation in acoustic travel times caused by 
the internal tides and internal waves is predicted. This will be eventually compared 
to actual perturbations in travel time recorded during the experiment. These pertur- 
bations are expected to be significant enough to be measurable. 

1.3 Thesis Content 

The scope of this thesis is as follows: Chapter two first presents the theory behind 
acoustic tomography, internal tides and internal waves. Next, using adiabatic normal 
mode theory, a formulation for the travel time perturbation of an acoustic signal due 
to internal waves is presented. Chapter three describes the region of the experiment, 
the experiment itself, and the data collected. Chapter four details the oceanographic 
data analysis conducted and reviews the assumptions which have been made. Chap- 
ter five examines travel time perturbations due to internal tides and internal waves. 
Comparisons between a simple ray model and the modal formulation derived in Chap- 
ter two are made. Finally, Chapter six presents conclusions and gives suggestions for 
future work. 
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Chapter 2 



Theory 



2.1 Acoustic Tomography 



Acoustic tomography can be broken into a forward and an inverse problem. In the 
forward problem, the travel time of a pulse along a path is computed from the sound 
speed and current fields of the ocean. In the inverse problem, given a set of travel 
times, the sound speed and current structure of the ocean is to be estimated. In this 
section, we will look at the ray theory approach to the forward problem. 

The formulation of the forward problem given here follows that of Howe et al. 
[1987]. Using geometric optics, the travel time T, along a ray path I\ in an ocean 
having a sound speed field C(x,<) and a current field u (x,f) is 



T - [ _ ( 

Jr, C(x, t) + u(x, t) ■ r’ 



(2.1) 



where s is the arc length, and r is a unit vector tangent to the ray. The sound speed 
field C(x, t) can be decomposed as 



C(x,<) = C 0 (x) + 6C(x,<), (2.2) 

where Co(x) is a “known” background sound speed field and 6C(x, t) is a sound speed 
perturbation such that 6C(x,t) <C Co(x). From these equations it can be seen that 
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a ray traveling with a current has a smaller travel time than a ray traveling against 
the current. 

The perturbation travel times for a reciprocal transmission between transceivers 
can be found by substituting (2.2) into (2.1) and linearizing. This gives 



ST + 



T = T? - To , ;= - J r 



8C(x, t) + u(x, t) ■ r 



ds , 



(2.3) 



and 



6T t r = T--T oi =-J r 



8C{x, t) — u(x, t ) • r 



ds, 



(2.4) 



where the + superscript refers to a ray propagating in the +x direction and the 
— superscript is a ray propagating in the — x direction. By taking the sum s,-, or 
the difference d, of (2.3) and (2.4), either u(x, t) ■ r, or 8C(x,t) may be eliminated, 
respectively giving 

28C(x, t) 



si = 8T? + 8T~ = - [ 

Jr, 



C o(x) 



■ds, 



(2.5) 



and 



4 -a? -sir— l ( 2 - 6 ) 

Howe et al [1987] showed that internal waves and large scale velocity shears do 
not significantly perturb the reference ray path so that for a two-way transmission the 
ray paths are reciprocal to first order and for a one-way transmission the ray paths for 
the reference ray and the perturbed ray are identical. In the BSPF experiment, only 
the one-way transmission case is examined. By noting that 6C(x,t) u(x,2) • r, 

(2.3) can be reduced to 



8T+ = T,+ - T 0l = 



r <5C(x,f) 

Jr, Cq(x) 



(2.7) 
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2.2 Internal Motions in the Ocean 



Internal motions in the ocean of interest in this thesis are internal tides and internal 
waves. This classification is a matter of the frequency of interest. Internal tides are 
forced at tidal frequencies Wtido and are either propagating or evanescent depending 
on whether u> t ide is greater than or less than the inertial frequency /, respectively. The 
propagating internal wave band is between the inertial frequency / and the highest 
buoyancy frequency N max • These internal motions are clearly seen in Fig. 2-1, which 
is an example of a temperature record that has been has been Fourier transformed 
into its frequency components (using methods to be discussed later in Chapter 4). 
Clearly shown in the figure are the semidiurnal internal tide along with its harmonics, 
and the internal wave spectrum at higher frequencies. N max at this depth is 215 cycles 
per day. 

2.2.1 Internal Tides 

Internal tides are generated by the surface or barotropic tides advecting water over 
the ocean bottom. Any change in the bottom topography causes the isopycnals to be 
displaced. This disturbance excites the internal modes in the water. The continental 
shelf break is a principal generation site of internal tides, whose generation mechanics 
have been well studied [e.g. Baines , 1974; Rattray et al 1969]. Internal tides can 
also be generated at other topographic changes such as guyots, troughs, sills, and sea 
mounts [e.g. Maxworthy , 1979]. 

Wunsch [1975] reviewed several internal tide records and noted they often showed 
a quasi-semidiurnal or quasi-diurnal periodicity, but that due to the record length, 
background noise, and frequency shifting of the tidal peak, it was often hard to resolve 
the tidal constituents. He also noted that the tidal signal was intermittent. Baines 
[1986] suggested that the intermittent nature of the internal tide might be due to 
the changes in stratification between the generation site and the observer, or due to 
frontal changes or upwelling. The generation of internal tides is a very non-linear 
effect. Many records show periodicities of 6 hours (twice semidiurnal) or 8 hours 
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Figure 2-1: Power Spectrum Density of Brancker temperature logger #3663 located 
on NE array at a depth of 32.1m. Evident are the M2 tide, its higher harmonics, and 
the internal wave band. N max at this depth is 215 cpd. The estimate has 4 degrees 
of freedom. 

(diurnal plus semidiurnal), or even 3 hours [ Baines , 1986]. The 3 hour (8 cpd), 6 
hour (4 cpd), and 12 hour (2 cpd) periods are evident in Fig. (2-1). Some of these 
higher harmonics of M2 could be due to computing the spectrum of temperature, 
rather than density, however, as will be discussed. 

In his review, Wunsch [1975] also noted that most of the internal tidal energy is 
contained in the first few lowest modes with mode one dominant. The direction of 
propagation is generally shoreward, but this need not be the case. For instance, Levine 
and Rickman [1989] performed a study which showed that the internal tide off the 
coast of California propagated along the shelf. They suggested that the internal tides 
were either being refracted from fronts along the coast or from mesoscale structures, 
or that the resultant direction is due to many waves being radiated from several shelf 
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breaks oriented in different directions. Another study by Holloway [1984] found that 
the first mode dominates and propagates onshore at an angle of about 30 0 from 
normal to the bathymetry. 

2.2.2 Internal Waves 

Internal waves are found throughout the ocean. They exist in the ocean interior 
due to disturbances between fluids of differing densities. Although numerous internal 
wave generation mechanisms have been detailed, only the more important mechanisms 
will be discussed here. There exist several thorough reviews which the reader may 
wish to consult [Thorpe, 1975; Muller and Olbers, 1975; LeBlond and Mysak , 1978; 
Olbers, 1983]. The internal wave generation mechanisms involve interactions with the 
bottom, the surface, or the interior. 

Bottom topography is an important mechanism in the formation of internal waves, 
particularly in coastal regions. This mechanism has been extensively studied in the 
region of the continental shelf break [Sawyer, 1983; Curtin and Mooers, 1975], over 
ridges in bays [Halpern, 1971; Chereskin, 1983] and over sills in straits [ Gargett , 1976; 
Farmer and Smith, 1980]. In a horizontally stratified ocean, the flow of the baroclinic 
tide over a sharp boundary (such as a continental shelf or a ridge) causes a shear 
instability and creates a lee wave. Depending upon the specific conditions, the lee 
wave will either propagate shoreward as a lee wave train or evolve into an undular 
bore structure. These internal wave packets or solitons have generally been observed 
when the tide turns. 

Surface interactions can also generate internal waves. Thorpe [1975] reviewed 
several surface interaction mechanisms, which we present here in order of assumed 
importance. First, the wind stress can produce a traveling stress field in the mixed 
layer, which results in a divergent current field in the mixed layer. The vertical 
Ekman velocities created by the wind stress can excite internal waves. Second, a 
traveling pressure field, created by atmospheric pressure, can resonantly couple with 
the internal waves to form waves that grow linearly. Third, a traveling buoyancy flux, 
produced by solar heating and precipitation, induces a vertical velocity at the lower 
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boundary of the surface layer. The induced velocities excite internal waves. Fourth, 
a spectrum of surface gravity waves can be produced by the wind stress acting over 
the ocean surface. Two of these higher frequency surface gravity waves, s 1 and s2 
can beat against each other to produce a lower frequency internal wave with wave 
number k,*, 

k^u; = k 5 i k 5 2> (2.8) 

and frequency u> tw 

Vi w = u 3 i - u s2 . (2.9) 

The interaction of the surface waves generates high frequency internal waves in the 
upper ocean, but this mechanism is insignificant for the deep ocean. Under conditions 
of a rough sea and a strong, shallow thermocline, however, it may be comparable to 
the wind stress mechanism (the first mechanism mentioned) [ Olbers , 1983]. Finally, 
instability between the surface boundary layer and the deep ocean, such as turbulence 
in the mixed layer, may also generate internal waves. 

As regards the role of the interior of the ocean in the generation of internal waves, 
very little is known. LeBlond and Mysak [1978] list as possible generation mechanisms 
processes such as penetrative convection, geostrophic adjustment, and internal tidal 
modulation of the thickness of the thermocline. 

Garrett and Munk [1972] looked at several internal wave observations and created 
from them a model of the internal wave field frequency-wavenumber spectrum. The 
Garrett-Munk model assumes horizontal isotropy and vertical symmetry. It has the 
form 

S{k h ,u>) oc /i _1 uT p+1 (u> 2 - / 2 ) -1/2 (2.10) 

where kh is the horizontal wavenumber, / is the local inertial frequency, and /i is the 
cut-off wavenumber which satisfies 

/i oc u/ -1 (u> 2 — / 2 ) 1/2 . (2.11) 

Substituting in (2.11) into (2.10), assuming u 2 / 2 , and then integrating (2.10) over 
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kh, gives for the frequency spectra 



S(w)=u;- p . (2.12) 

The indices p and r are determined from observed spectra, which suggested values 
of p between 5/3 and 2, and r between 2/3 and 3/2. Garrett and Munk proposed a 
value of 2 for both r and p. The Garrett-Munk model has been modified and refined 
over the years as new experimental information has become available. 

While the spectrum for the deep ocean is well modeled by Garrett-Munk, that 
for shallow waters is not. An internal wave experiment, conducted as part of the 
GARP Atlantic Tropical experiment or GATE in 1974, gave a spectrum which devi- 
ated strongly from the Garrett-Munk model [ Kdse and Siedler, 1980]. The frequency 
spectrum showed three dominant energy bands centered about the near-inertial fre- 
quency, the M2 tidal frequency, and 3 cycles per hour (cph). The high frequency 
wave field was intermittent and anisotropic. The non-stationarity of the spectra is a 
result of the variability of the atmospheric generation mechanisms [Kdse and Clarke , 
1978]. Kase and Clarke also argue that the the higher frequency waves are trapped 
near the surface where a maximum N overlies a thick weakly stratified lower layer. 
These waves form a low order standing mode. Baines [1986] also noted a study in 
which for depths less 800m, 25-45% of the energy was in mode 1 which propagated 
shoreward. 

Dynamical Equations for Internal Waves 

Using the basic equations of momentum and continuity, an exact modal solution for 
internal waves can be derived. The equation of momentum for an isentropic fluid in 
the absence of viscous effects is given by [e.g. Gill, 1982] 

Du , 

— + 2Dxu=-p- , Vp-g. (2.13) 
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and the continuity equation for an incompressible fluid is given by 



V • u = dujdx -f dv/dy -F dw/dz = 0. 



(2.14) 



Following the derivation of Apel [1987], a small perturbation in p and p is applied to 
(2.13). The perturbation pressure and density are then given respectively by 



P = -gpz + P, 



(2.15) 



and by 



P — po + p' ■ 



(2.16) 



The momentum equation (2.13) is linearized by neglecting the products of perturba- 
tion quantities, and becomes 



du 1 dp' 

dt V po dx' 


(2.17) 


dv 1 dp' 


(2.18) 


dw 1 dp' p' 

dt p 0 dz po 9 ' 


(2.19) 


where / is the inertial or “local Coriolis” frequency. We now assume a separable 


solution for the vertical velocity w(x) of the form 




w(x, t ) = w(z) exp[t(£\r + ly — )], 


(2.20) 


where w(z) is the vertical eigenfunction. The vertical velocity 
boundary conditions of no normal flow at the bottom 


tu(x) satisfies the 


te(x) = 0 at z = h (bottom), 


(2.21) 



21 



and that a particle on the free surface remains on it 



w(x) = drj/dt at z = 1/ (surface), 



( 2 . 22 ) 



which can be applied at z — 0 since 77 is small. This approximation introduces little 
error, as the vertical disturbance at the air-water boundary is much smaller than the 
maximum amplitude by an order of 1000. By manipulating equations (2.14-2.19), 
one can obtain 



d 2 w 



dz 2 



+ kl 



N 2 (z) — u> 2 



LO 2 -} 2 



w = 0, 



(2.23) 



where kh is the horizontal wavenumber that satisfies 



kl = k 2 + / 2 , 



(2.24) 



and where N is the buoyancy or Brunt-Vaisala frequency which satisfies 



N 2 = 



9 dp Q 
Po dz 



(2.25) 



For u> < f or u > N, the vertical wave number is imaginary and the solutions to 
(2.20) are evanescent. 

By using the relationship 



w — 



d A 

dt' 



where C is the particle displacement, and substituting in (2.23) gives 



dz 2 



+ k 2 



N 2 {z)-u 2 



LO 



2 _ 



P 



C = 0, 



with the following boundary conditions: 



(2.26) 



(2.27) 



£(2) = 0 at 2 = 0 (surface); 



(2.28) 



and 



£(2) = 0 at 2 = h (bottom). 



(2.29) 
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This equation can be formulated as an eigenvalue equation and is easily solved using 
numerical methods such as finite element methods. The exact modal solution is 

oo 

C = A ‘M Z ) ex P b (h x + by - w<)] , (2.30) 

;=i 

where l = 1,2,... is the mode number, Ai is the amplitude, and d>i(z) is the vertical 
internal wave mode eigenfunction. 



2.3 Adiabatic Acoustic Normal Modes 

We now return to acoustics. In a range independent environment, the acoustic pres- 
sure p(z, x) for a point source at r = 0, 2 = z\ is described by the Helmholtz equation 

V 2 p + I< 2 (z)p = -~6{z - Z\ )6(r). (2.31) 

r 

Brekhovskikh and Lysanov [1991] give as the solution 



p(r,z) = 7n <f > m(z:)<f>m(zi)HQ 1 \K m r). 

m 

The vertical mode function <p m satisfies 



(P (f>m 
dz 2 



+ [A' J (.-) - A'i] 4,„ 



= 0, 



(2.32) 



(2.33) 



where K(z ) is the total wavenumber and is the m th modal eigenvalue. Equa- 
tion (2.32) is a good approximation for many parts of the ocean, but as it is range 
independent, it is inadequate where the depth of the ocean changes significantly, 
where the sound waves cross ocean fronts, or where there is a significant change in 
the sound speed profile along the propagation path. In these cases the Helmholtz 
equation, equation (2.31), must be modified to take into account the variation of 
ocean parameters with range and becomes 



V 2 p + I\ 2 (z, r )p = 0. 



(2.34) 
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Shang [1989] gives as the solution to (2.34) 



p(r,z) 



E 



F n {r)4> n (z,r ) 



(2.35) 



where F n (r) is the radial component and the vertical mode function <f> n (z,r) satisfies 



8 ( 1 dtn\ ( K*{z,r) K(r)\ 

dz \p 0 (z) dzj^y po(z) po(z) ) 9n • 



(2.36) 



The mode functions <f> n are mutually orthogonal and satisfy 




(2.37) 



Substituting (2.35) into (2.34), multiplying by <f> m , and then integrating over 2 from 
0 to h results in the following coupled equation 

+ ('CM + ^) F„ = -Y, [A™(r)F n + 2B m „(r)J (2.38) 



where A mn and B mn are coupling coefficients defined as 



A mn = J p 0 1 <f> m {r,z)^j-(z,r)dz, (2.39) 

B mn = J po 1 <f> m (r,z)^-(z,r)dz. (2.40) 

To obtain the “adiabatic approximation”, we neglect the l/4r 2 term in (2.38) since 
we are interested in solutions with K m r 1, and we assume that the mode coupling 

effects are small. Physically, this means that the ocean properties (depth, sound 

speed) “vary slowly” as a function of horizontal range. The phrase “vary slowly” 
means that the fractional change of the ocean properties over a distance of less than 
a wavelength is small [ Pierce , 1965]. This type of medium has been termed “almost 
stratified”. The adiabatic approximation solution to (2.34) is 



p[r,z) = V2xexp(i7r/4)2^ . y 

rn yj K m (r )r 



exp 




(2.41) 
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In the Barents Sea, the background sound speed profile and the depth vary rea- 
sonably slowly along the sound propagation path on either side of the front. Taking 
advantage of this we can crudely treat the path on either side of the front as a “nearly 
range independent” segment, and use the adiabatic approximation. We expect strong 
coupling of the modes, though, in going through the frontal region. 



2.4 Sound speed perturbations due to Internal 
Waves 

Equation (2.3) shows that the two mechanisms affecting the perturbation of the acous- 
tic signal are the horizontal current u(x, t) ■ r and the sound speed perturbation 
8C(x,t) which is related to the particle displacement £ using the following relation- 
ship 

sc = ^(, (2.42) 

where C' = dC v /dz is the potential sound speed gradient and ( is the solution to 
equation (2.27). Assuming that both mechanisms are due to internal waves, Essen et 
al [1983] showed that SC has a much greater affect on the acoustic signal than does 
u(x, t) • r. 

Following a recent derivation by Jin [1992], for a band of frequency a;, equa- 
tion (2.30) can be generalized to 

00 [N r2ir 

C = YL / «/(^, 0)t/>/(u;, z) exp[i(ki(u;)x + U{u)y — ut)\dwd9, (2.43) 

where ai(u,0) is the internal wave amplitude spectrum for mode /, il>i(u,z) is the 
internal wave mode function, 0 is the azimuthal angle, and ki(u>) and Ii{uj) are the in- 
ternal wave eigenvalues that satisfy equation (2.24). The quantities il>i(u>,z), ai(u>,0), 
ki(u), and //(a;) are assumed to be range independent. 

We next follow Rajan et al. [1987] and perturb the background acoustic wavenum- 
ber I\ °(r, z) by perturbing the sound speed profile. A new wavenumber K(r,z) can 
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then be written as 



K (r, -| = w/[C 0 (r, z) + 5C{r,z)}. 



(2.44) 



Using standard perturbation theory, a perturbation eigenvalue SI\ m (r) for (2.41) is 
found to be 



roo 

6K m {r) = (Q) I 4>W{r,z)\&K{r,zf\4> ( °\r,z)dz. (2.45) 

/ A li ( r J 0 



2 Km (r) 

The middle factor in the integral in (2.45) can be rewritten as 



81<{r,z) 2 = I< 2 {r,z)~ h'W 2 {r,z) 



u 



C2(r, z) 



-2 



' fCfrvz) 

„ Co(r,z)J 



'= -K m7 (r,z)(26C/C 0 ) 



(2.46) 



Substituting (2.46) into (2.45) results in 



6Km(r) = dfc jf'" w i^*>r < 2 - 4? > 



Using the WKB approximation, the adiabatic modal phase is 



,R 

<p{r) = / I< m (r)dr, 
Jo 



(2.48) 



and the background modal phase is 



¥>V) = I" Kl(r)dr = K° n , 

Jo 



(2.49) 



The modal phase perturbation is thus 



8<~p m = / 8I\ m (r)dr 
Jo 



(2.50) 
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If the received signal is processed by a narrow band filter with center frequency 
u> 0 , then the group travel time of each modal pulse is given by the stationary time as 



duo 



so that the group travel time perturbation is 



dtm — 



d[d(y9 m ] ufp 

du 



(2.51) 



(2.52) 



Inserting (2.47) into (2.50), and (2.50) into (2.52) and using the results of equa- 
tion (18) of Lynch et al. [1991] results in 



St m = / / G m (r,z)6C{r,z)drdz, 

Jo Jo 



(2.53) 



where G' m (r, z), the group velocity dispersion, is 



G m {r, z) — pQ 





Ct 0| (r) 

u£\t) 








, (2.54) 



and C£\r) and U£*(r) are the modal phase and group velocities, respectively, for 
the unperturbed case. By assuming range independence, and also that 




(2.55) 



(2.54) becomes 

\t£\z)\ 2 + u 0 ^\<j)£ ) {z )\ 2 • 

where p(z ) = 1 and ~ Co- The bar denotes the average sound speed. 
The covariance of 6t m is formed by 



Gm{z) ^2 



(2.56) 



cov m (R,T ) = (6t m (R,t)8t* m (R,t + r)), (2.57) 
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which using (2.53) becomes 



rR rR r oo r oo 

cov m {R,r) - / / / / G m (r, 2 )G m (r / ,z') x 

JO Jo JO JO 

(6C(r, z, t)6C(r\ z', t + T))drdr dzdz ' , (2.58) 

where from (2.42) and (2.43) 

00 /-A/ ,-2>r 

6C(x,y,z,t ) = C' p (z)J2 / / <*/(<*>, 0)0/(u>, *) x 

i J f Jo 

exp[i(^;(u;)a: cos(0) + khi{oj)y sin(0) — ut]du>d9. (2.59) 

Assuming that the amplitudes ai(oj,0) are independent, random variables with zero 
means implies that 

(a,(u,9)) = 0, 

and 

(ai(u,0)c if,(w,tf)) = | a t (u,0)\ 2 S^'dee'Sir, 

where 8 WW ’, Sg$>, and 6u> are three delta functions. Additionally, if it is assumed that 
the x axis is in the direction from the source to receiver, we obtain 

(■ 6C(r,z,t)6C(r\z’,t + T)) = 

00 r N r 2n 

Y / \a,(u,0)\ ip,(u,z)tpj(u',z') x 

1 = 1 ** f 

C p (z)C p (z') exp [i(A;w(a;) cos 0(r — r ') — wr)] dudO. (2.60) 

Substituting (2.60) into (2.58) results in 

00 t N r2n 

cov m (R, T ) = J2 / / \ai(uj,0)\ 2 pi{u,0)e~' u ’ T dudT, (2.61) 

i=i J f 
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where 



r oo r oo 

p,(u,0) = 11 M^^:(u\z f )C' p (z)C' p (z')x 

\lo Jo z )G m (r', Z ') exp[z^/(w) cos 0(r - r')]drdr' 



(2.62) 



dzdz . 



If G m { r ,z) is assumed to be range independent, then it can be taken outside the 
integral and (2.62) reduces to 



rR rR r oo 

p,(u,0)= / / exp[tfc w (u>)cos0(r - r'))drdr' / ^,(w, z)C p (z)G m {z)dz 

Jo Jo Jo 



(2.63) 



The first integral on the right hand side of (2.63) reduces to 

2[1 — cos (Icm cos OR)] 



and can be further reduced to 



(khi cos #) 2 



r>2 • 2 ( ^hlR J\ 

R sine I — — cos v , 



where sinc(x) = sinx/x. Thus (2.63) becomes 



Pi(lj,0) = R? sine 2 ( — cos 9 



r oo 

J q ipt{u,z)C' p (z)G m (z)dz 



(2.64) 



The covariance spectrum can be represented by 



/ OO 

cov m (R , T)e~' UT dr. 

-OO 



(2.65) 



Inserting (2.61) into (2.65) yields 



F m (R,0) = J Q |ai(^,0)| 2 p/(w,0) (y e ,(n u ’ )T dr) dudd, (2.66) 
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which after integrating over r gives 



00 rN r 2 jt 

F m (R,n) = '£ / \a,(^d)\ 2 Pl (uj,e)6(n-u)(L J d0, 

l— 1 J f ’'O 

which can be reduced to 

00 r 2 * 

F m (R,n)=J2 / \ai(Sl, 0 )\*p,(Sl y O)dO. 

1=1 J ° 

The previous results can be simplified by assuming that 

| a,(u> y 6)\ 2 = A(u)B(9)D h 



and further simplified by assuming B(9) = 6(6 — 0 O ) and that the packet 
internal waves is Rd . For the case of 9q = 7r/2, (2.61) becomes 



cov r 



i(i?)r=0 — Rd f 

J J 



ZMWi 



[ H MV, z)C' p (z)G m (z)dz 

Jo 



dCl. 



In this thesis we will evaluate equations (2.53) and (2.70). The assumptions 
in evaluating those equations are detailed in Section 4.3. 



(2.67) 



(2.68) 



(2.69) 
width of 



(2.70) 
we make 
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Chapter 3 



Barents Sea Polar Front 
Experiment 



The Barents Sea Polar Front (BSPF) experiment took place from August 6 to August 
26, 1992 in the coastal waters off of Svalbard, about 42 nautical miles due east of 
Bjprnpya or Bear Island at the northern edge of the Bj0rn0ya Trough. The BSPF 
experiment involved three CTD and ADCP surveys, and the deployment of four 
moorings for about ten days. Water depth in the experimental area varied from 
1 10 m to 380 m. 



3.1 Barents Sea 

The Barents Sea is a shallow, peripheral basin bounded by Svalbard on the north, 
Norway and Russia on the south, Novaya Zemlya on the east, and the Norwegian Sea 
on the west. The majority of the sea is contained between 70° and 80° north latitude 
and 16° and 60° east longitude (Fig. 3-1). The boundary between the Norwegian 
Sea and the Barents Sea is the transect from Spitsbergen to Bear Island to the coast 
of Norway. Table 3.1 gives the characteristics of the waters discussed below. For a 
detailed description of the water masses and types see Hopkins [1991] or Johannessen 
[1986]. The Barents Sea Polar Front (BSPF) is formed at the boundary of two 
currents, the North Cape Current to the south, and the Bear Island Current to the 
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Figure 3-1: Chart of the Barents Sea with major currents. The thick, solid line is the 
Barents Sea Polar Front. 
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Table 3.1: Water Mass Definitions 



Water Mass 


Symbol 


T(°C) 


S'(ppt) 


Comments 


North Atlantic Water 


NAtW 


> 8 


> 35.3 


Surface water of 
N. Atlantic Current 


Norwegian Atlantic Water 


NwAtW 


> 2 


> 35.0 


Water of Norwegian 
Atlantic Current 


Polar Water 


PW 


Freezing 
to 2 


33.0-34.4 




Barents Polar Water 


BrPW 


-1-3 


< 34.3 





north. The North Cape Current, an extension of the Norwegian Atlantic Current 
(NwAtC), carries warm, saline Norwegian Atlantic Water (NwAtW) into the Barents 
Sea. The Bear Island Current, which flows along the northern part of the Bjprnpya 
or Bear Island Trough, carries cold, fresh Barents Polar Water (BrPW) out into the 
Norwegian Sea. The front is a stable interface because the warm, saline NwAtW 
has approximately the same density as the cold, fresh BrPW. The BSPF is some 
1500 kilometers long and roughly follows the 200 meter isobath. It is also strongly 
influenced by the tidal cycle and can be displaced by a much as 10 kilometers during 
a tidal cycle [Johannessen 1986]. 

The tides in the Barents Sea are primarily semidiurnal (M 2 ) with a weaker diurnal 
component (Ki). Both Gjevik [1990] and Schwiderski [1986] have proposed tidal 
models for the Barents Sea. The predicted amplitudes for each model are given in 
table 3.2. The models differ in that Gjevik’s model resolves amphidromes for both the 
M 2 and K\ in the region of Bear Island, whereas Schwiderski’s does not. Huthnance 
[1981] made a series of current meter measurements from Northern Norway to Bear 
Island. From his experiment, he estimates that both the M 2 and K\ are mostly 
barotropic (about 86 % ). The predominantly barotropic S 2 tide is smaller than the 
M 2 , and has a fairly consistent relation with it. 
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Table 3.2: Tidal Amplitudes at Bear Island 





Obs. 


Gjev 


Schw 


Constituent 


Amp. 


Amp. 


Amp. 




(cm) 


(cm) 


(cm) 


a/ 2 


34.2 


34.5 


40 


Ki 


5.4 


6.3 


5.5 



Obs. — Observed values taken from Gjevik [1990] . 
Gjev — Values taken from Gjevik [1990] 

Schw — Values taken from Schwiderski [1990] 



3.2 Experimental Objectives 

The purpose of the experiment was to study the structure and the dynamics of the 
Barents Sea Polar Front by using a variety of moored and shipboard oceanographic 
instruments. A key aspect of the BSPF experiment was the Barents Sea Tomography 
Transmission Test (BSTTT) whose major goals were: first, to determine the iden- 
tifiability, resolvability, strength and variability of the multipath structure between 
the acoustic moorings; second, to examine the acoustic signature of various oceano- 
graphic processes (surface waves, internal tides, internal waves, eddies, etc.) between 
the acoustic moorings; and third, to use “tomographic inverse images” to study the 
region of the BSPF. 



3.3 Data Acquisition 

Four acoustic instruments were deployed during the experiment in a trapezoidal shape 
about 50 kilometers apart at the base, with the distance between the northeast and 
southeast moorings being 34 kilometers (Fig. 3-2). The Northwest (NW) and South- 
west (SW) moorings contained 400 Hz tomography transceivers while the Northeast 
(NE) mooring contained a 224 Hz projector. These three moorings were deployed for 
a period of 9 days (NE) to 12 days (SW). The transceiver on the NW mooring acted 
as 400Hz receiver only, while the transceiver on the SW mooring acted intermittently 
as a 400Hz source and continuously as a 400Hz receiver. The fourth or Southeast 
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74.75 




Figure 3-2: Chart of the BSPF experiment. Triangles are array positions. Arrows 
denote the acoustic paths during the experiment. 

(SE) mooring had a vertical line array(VLA) of sixteen hydrophones which collected 
70 hours of data. The data from the VLA was telemetered in real time back to 
the research ship using a wireless local area network. The acoustic paths between 
the moorings are shown in Fig. 3-2. Details of the vertical line array are given in 
von der Heydt et al. [1992]. All of the moorings except the SW mooring contained 
several current meters and Brancker temperature loggers to provide a vertical profile 
of current and temperature at each mooring. Details are given in Table 3.3. 

Mooring motion of the source or receiver due to tides or other currents changes 
the source-receiver range, and introduces error into the measurement of the travel 
time. To correct for this mooring motion, all of the moorings except for the NE 
mooring used acoustic navigation systems. The determination of mooring motion is 
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Figure 3-3: CTD grid employed during BSPF experiment. Circles are CTD stations. 
The grid was fully covered three times. 

accurate to within lm. Additionally, the SE mooring had two tilt meters. The ship 
also conducted three complete 70 km by 80 km CTD grids (Fig. 3-3) of the research 
area, two detailed grids across the Barents sea polar front, and a 24 hr CTD anchor 
station at the SW mooring. Details of the CTD grids and the anchor station are given 
in Pawlowicz [1992]. Several SUS charges were also dropped to allow determination 
of bottom properties and reverberation characteristics. 
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Table 3.3: Barents Sea Moorings 



Northwest Mooring-170m 



Northeast Mooring-142m 



Southwest Mooring— 380m 



Southeast Mooring— 280m 



Instrument 


Depth (m) 


Neil Brown SACM 


19.4 


Brancker Tlog 


29.5 


Brancker Tlog 


39.2 


Aanderaa CM 


49.3 


400 Hz transc. 


66.7 


Aanderaa CM 


159.9 



Neil Brown ACM 


22.1 


Brancker Tlog 


32.1 


Brancker Tlog 


42.1 


Neil Brown ACM 


52.0 


Aanderaa CM 


82.1 


224 Hz proj. 


121.8 



Neil Brown SACM 


19.9 


Brancker Tlog 


29.9 


Brancker Tlog 


39.9 


Aanderaa CM 


50.0 


Aanderaa CM 


80.0 


Brancker Tlog 


140.0 


Brancker Tlog 


200.5 


400 Hz transc. 


358.1 


Vertical Array 


126-275.8 
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Chapter 4 



Temperature Field Data Analysis 



The raw data from which the calculations in this chapter are made comes mainly from 
the CTD casts made at each mooring site, and the moored Brancker temperature 
loggers. Some temperature information is also from sensors in the current meters. 
Figure 4-1 is a typical Brancker temperature logger record. The temperature record 
is non-linear in that at higher temperatures the signal is “spikier” than at lower 
temperatures. This effect is most pronounced for those loggers located where the 
temperature gradient is most rapidly changing. A one meter downward displacement 
of water at the “knee” of the temperature curve (i.e. the mixed layer base) results in 
a much greater change of temperature than a one meter upward displacement. Figure 
4-2(a) shows a sinusoidal 9.9m amplitude particle displacement which, acting on a 
typical temperature profile (Fig. 4-2(b)) results in a “distorted” temperature record 
(Fig. 4-2(c)). 

The changes in temperature seen are due to both a vertical particle displacement 
of water caused by internal tides and waves, and a horizontal advection of water 
caused by the barotropic tide. Which mechanism is dominant can be determined 
by evaluating constraints imposed by the local temperature gradient, current veloc- 
ities, and water depth. Figure 4-3 shows three CTD temperature profiles taken at 
the Southwest array. If we take profile #418 as the baseline, and translate profile 
#218 vertically by +7m, then profile #218 will lie directly over profile #418. This 
corresponds to a temperature change solely due to internal tides and internal waves. 
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Figure 4-1: Temperature record from Brancker logger # 3666 located on the SW 
mooring at a depth of 29.9m. Brancker temperature resolution is 0.1 °C 

Similarly, if profile #118 is translated upwards in depth by 15m and moved hori- 
zontally by 0.25 °C it will lie over profile #418. This corresponds to a temperature 
change due to both a barotropic current and internal tides and internal waves. We will 
now consider the different effects separately, and thereby determine the percentage 
that each signal contributes to the whole. 

4.1 Horizontal Advection versus Internal Tides 
and Waves 

The amount of horizontal advection was determined at each of the three instrumented 
arrays by using the assumption that the temperature changes were solely due to 
horizontal advection and that any vertical motion due to internal tides and waves 
or topographic changes was negligible. In order to explain the method, the data 
from the NE mooring will be used. Figure 4-4 shows the temperature and velocity 
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Figure 4-2: (a) is a 9.9m sinusoidal particle displacement about the logger depth 
(dashed line), (b) is a smoothed temperature profile at the SW mooring. The dashed 
line is the logger depth, (c) is the “distorted” temperature record which results from 
the sinusoidal particle displacement in a varying temperature gradient. The dashed 
line is the axis of zero phase. 
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Figure 4-3: Three temperature profiles at the SW array taken over a span of 18 
days. The temperature records show both a temperature change due to vertical 
displacement (caused by internal tides and waves) and due to horizontal advection 
(caused by the barotropic tide) 

records from Aanderaa #9259 at a depth of 82.1 meters. The temperature record has 
a minimum of —0.9 °C and a maximum of —0.36 °C with a peak to peak change in 
temperature of 0.208 °C. 

Figure 4-5 is a contour map of temperatures taken from the CTD profiles sur- 
rounding the NE mooring and shows that the warmer water at 82 meters is to the 
east and to the south, with the east-west gradient being much larger than the north- 
south gradient. Several things are apparent in Fig. 4-4. The u velocity lags the 
temperature record by 90 degrees, while the v velocity is in phase with the tempera- 
ture record. Given that the u velocity is over twice the v velocity and given that the 
east to west temperature gradient is much larger than the north south temperature 
gradient, the predominant effect will be from the u velocity. Since the u velocity lags 
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Figure 4-4: U and V velocity, and Temperature records from Aanderaa #9259 at the 
NE array at a depth of 82.1 meters. 




Figure 4-5: Temperature contour map at NE mooring at a depth of 82.1m using 
adjacent CTD stations. The current ellipse from Aanderaa # 9259 shows the amount 
of horizontal advection during one tidal cycle. 
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temperature by 90 degrees, at the time when the current shifts from flowing westward 
to flowing eastward, the temperature will be at a maximum, i.e. it is at this time that 
the warm water to the east will have been advected the greatest distance to the west. 
This can also be seen by taking the integral of the u velocity; the displacement record 
is now 180 degrees out of phase with the temperature record which is consistent with 
having the warmer water to the east. 

The displacement ellipse was then plotted over the temperature contour. From this 
plot, one sees that the maximum change in temperature due to horizontal advection 
is 0.037 °C, which only accounts for 18 percent of the peak to peak temperature 
change. Another argument against horizontal advection being the mechanism behind 
the temperature changes is that the —0.9 °C water would have had to been advected 
from over 14 km away from the northwest and the —0.36 °C water would have had to 
been advected from over 10 km away from the east or over 20 km away from the south. 
None of these distances are credible given the results listed in Table 4.1. Analyses of 

Table 4.1: Horizontal Advection of Water due to Barotropic Tides 



Array 


Depth (m) 


N-S disp (km) 


E-W disp (km) 


NE 


82.1 


0.98 


2.2 


NW 


49.3 


1.76 


2.46 


SW 


80.0 


0.89 


1.82 



the SW and NW moorings yield similar results. 

4.2 The Internal Tide and Internal Wave Spec- 
trum 

In this section we convert temperature fluctuations into particle displacement fluctua- 
tions and then into internal tide and internal wave amplitudes using Fourier analysis. 
In performing this analysis we assumed that the temperature fluctuations can be 
used to obtain the magnitude of the particle displacement. This will be true if the 
temperature-salinity (T-S) relationship is constant, i.e. one can go uniquely from 
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Figure 4-6: Temperature and Salinity versus depth for the CTD casts taken at the 
NE array. 

temperature displacement to density displacement. Figure 4-6 shows the relative 
constancy of the T-S relationship over three different CTD casts at one site. A rise 
in temperature for a sensor (see Fig. 4-7) above the mixed layer (i.e. negative tem- 
perature gradient) corresponds to a downward displacement of warmer water. The 
vertical displacement Sz for a parcel of water is 

^ = ( 4 -l) 

where 8T is the temperature fluctuation measured by the Brancker logger, and dz/dT 
is the temperature gradient at the mooring. We have taken z as being negative 
downward. The particle displacement can be determined by integrating (4.1) so that 
each temperature corresponds to a given depth. This method fails if the temperature 
profile is double valued near the depth of the logger or the logger is near the knee 
of the T(z) profile. Using the temperature profiles from the CTD casts, we obtained 
a best fit temperature profile which preserved the slopes of the CTD temperature 
profiles. Each temperature record was then converted to a particle displacement 
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Figure 4-7: Temperature record for Brancker logger #3663 at NE mooring at a depth 
of 32.1m. 



record (see Fig. 4-8). The method assumed that the best-fit linearized temperature 
profile was temporally stationary, whereas in fact it is not quite, as demonstrated 
by Fig. 4-3. Additionally, in several records the temperature minima and maxima 
recorded by the Brancker temperature loggers exceeded the temperature minima and 
maxima recorded during the CTD casts. The CTD casts represent “snapshots” of the 
temperature profile, and of the ten casts taken at the NE, SW, or NW moorings only 
two occurred (one at the NW mooring and one at the NE mooring) when the moorings 
were deployed, and these CTD casts were at the ends of the temperature records. In 
order to compensate for this mismatch, the temperature profiles were extended at the 
upper and lower ends while maintaining the same slope. This process was analogous 
to including some horizontal displacement. It must be remembered that there are 
two processes taking place, horizontal advection and vertical displacement, and that 
when the Brancker temperature exceeds the CTD temperature, what has occurred is 
obviously both a horizontal advection and a vertical temperature displacement. 

To estimate the error of extending the temperature profiles, the temperature 
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Figure 4-8: Particle displacement £ calculated for Brancker logger # 3663 on NE 
array at a depth of 32.1m. 



record for Brancker #3660 was examined. This temperature record, like the others at 
the SW mooring, showed several eddies, which act similar to a horizontal advection of 
the profile (e.g. Fig. (4-1). The particle displacement and amplitudes were computed 
two ways for the Brancker #3660 temperature record: first, assuming that the eddy 
was a horizontal advection and adjusting the temperature profile accordingly; and 
second, filtering out the eddy feature and using an unadjusted temperature profile. 
The difference between the amplitudes in the internal tide band computed the two 
different ways was at most less than 3.5%. 

Once the particle displacement was computed, it was then Hanning windowed, and 
analyzed using the Welch method of fast Fourier transform(FFT) spectrum analysis. 
Figure 4-9 is a typical spectrum. The most significant peak in the internal tide region 
is the semidiurnal peak. The amplitude was extracted by summing the area under 
the peak 
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Figure 4-9: Power spectrum density for the particle displacement calculated for 

Brancker logger # 3663 on the NE array at a depth of 32.1 m. The estimate has 
four degrees of freedom. 

where A is the amplitude, PD(k) is the power density function at k and N is the length 
of the sequence. Table 4.2 gives the amplitudes at the Southwest array computed 
using (4.2) for the M2 tidal component. 

Due to the shortness of the records, high resolution Fourier analysis is not possible. 
This results in combined tidal amplitudes for some frequency values. The M2 peak 
probably contains S 2 and other semidiurnal constituents, though they are small. Also, 
the amplitudes computed this way are average amplitudes, and therefore smooth out 



Table 4.2: SW Mooring Particle Displacement Amplitudes £ 



Logger Depth(m) 


M2 (m) 


29.9 


1.8 


39.9 


2.1 


140.0 


7.9 


200.5 


6.1 
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Figure 4-10: Amplitude spectrum for Brancker logger # 3660 on SW array at a depth 
of 39.9m. 

the previously noted intermittent nature of the internal tides and internal waves. 
This method gives reliable estimates for frequencies in the internal tide range (1-10 
cpd), but may underestimate the amplitudes in the higher frequency internal wave 
band. This effect will be addressed when estimating the amplitude spectrum used in 
equation (2.70). Figure 4-10 shows a typical amplitude spectrum. 

The amplitudes extracted using Fourier analysis are used to determine the relative 
weights of each of the normalized internal tide and internal wave modes. This analysis 
is performed for the SW mooring because it contains four Brancker temperature 
loggers and spans 170m of the water column, while the NE and NW moorings only 
have two temperature loggers and span only 10m of the water column. Also, using 
a least squares fit, up to four modes can be estimated at the Southwest array (Fig. 
4-ll(a)). The least squares fit was simplified by averaging together the two shallow 
loggers, as they are only 10m apart and can be approximated as a single sample. 
Figure 4-1 1(b) is the result of this average. The results from both methods were 
consistent. The fit was only performed for frequencies up to 30 cpd, since this is 
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the maximum buoyancy frequency for the two deep loggers, and in fact, the method 
appears to break down at about 27cpd. Also, the sign of the mode shapes was 
consistent in the internal tide band (0-8 cpd), but the sign of the third mode changed 
in the frequency band of 10-30 cpd. This change coincided with change from the 
logger at 140m having the largest amplitude in the internal tide bank to the logger at 
200m having the largest amplitude in the internal wave band. The percentages were 
averaged over the internal tide and internal waves, resulting in an average percentage 
for each of the first three modes of: mode #1 - 0.58, mode #2 - 0.26, and mode #3 
- 0.16. As a check on these numbers, the modes were also estimated incoherently, by 
performing a least squares fit of the first mode squared to the square of the amplitudes. 
The residual was then least square fitted to the second mode squared. This process 
was performed for the first three modes at the M2 tidal frequency. The resulting 
relative percentages for the first three modes were: mode #1 - 0.62, mode #2 - 0.22, 
and mode #3 - 0.16. Both of these sets of numbers are consistent with each other 
and the results reviewed in Wunsch [1975]. 



4.3 Review of Assumptions 

We are now ready to compute the acoustic fluctuations using Eqs. (2.53) and (2.70), 
but before doing so a review of the validity of the assumptions is warranted. First, 
we have assumed, in truncating the IW mode series to three modes, that the internal 
tide and internal wave spectra are red i.e. most of the energy is in the lowest modes. 
This is probably an excellent assumption, as numerous experiments have concluded 
that the bulk of the energy in both the internal tide and internal wave bands is in the 
lowest modes. The ratios of the first three internal wave modes calculated in Section 
4.2 show good agreement with those in Wunsch [1975]. 

Second, we have assumed in solving equation (2.70) that the direction of propaga- 
tion is along-shelf instead of across-shelf. We assumed this direction of propagation 
for two reasons: first, we didn’t have available with our single point moorings any 
information on the direction of propagation; and second, this direction ( 0 q = 7r/2) 
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Figure 4-11: (a) is the least squares fit of four amplitudes measured at the SW array 
to a sum of the first four IW modes, (b) is the same calculation except that the 
two shallow Branckers are averaged together. In both figures, considerable noise is 
evident in the internal tide region, and the method is unreliable above about 27cpd 
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maximizes the sine 2 function and thereby, the travel time perturbation covariance 
(Eq. 2.61). This gives the largest effect we could hope to see, and thus represents 
an “upper limit” calculation. Additionally, although the assumption of along-shelf 
propagation is more tenuous than that of across-shelf propagation, it is perhaps not 
wholly unjustified, since the experimental region is on the edge of the Bear Island 
trough and in the vicinity of the strong Barents Sea Polar front. As discussed in 
Section 2.2.2, although internal tides and internal waves propagate shoreward, one 
experiment noted that the internal tides propagated along the coast. It was postu- 
lated that either coastal fronts, or mesoscale structures were refracting the internal 
tide or that internal tides radiating from differently oriented shelves were causing a 
complicated wave structure [Levine and Rickman, 1989). The experimental area con- 
tains both such features. The Bear Island trough certainly qualifies as a shelf area 
with different orientations. 

Third, we have assumed that the change in temperature seen by the logger was 
largely due to a vertical particle displacement caused by internal tides and internal 
waves rather than due to a horizontal advection caused by the barotropic tide. We 
showed in Section 4.1 that at most 18% of the change in temperature could be at- 
tributed to horizontal advection. This effect was included implicitly in determining 
the particle displacements, but by not including it resulted in a less than 3.5% error. 

Fourth, we have assumed that 8C(x,t) u(x, t) • r. From the Aanderaa data 
the maximum current was 0.325 m/s, whereas in the mixed layer, a 5m particle 
displacement could result in a 6.0 m/s change in sound speed. 
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Chapter 5 



Travel Time Disturbance 
Predictions 



Using the theory from chapter two, and the internal wave frequency spectrum and 
amplitudes generated using the methods of chapter four, a modal picture and a ray 
picture of the travel time perturbations are now generated. These predictions will 
eventually be compared with the actual results from the experiment, when these 
become available from colleagues of the Naval Postgraduate School who are working 
on the travel time estimation problem. 

5.1 Theoretical Results for the M2 Tide 

In chapter two, a formulation for the travel time perturbation, equation (2.53), and 
travel time covariance, equation (2.70), was derived, using adiabatic normal modes as 
a basis. The adiabatic assumption requires a “slowly” varying sound speed profile and 
bathymetry. We further assumed range independence for the group velocity dispersion 
G m (z ) and for the particle displacement (. In evaluating equations (2.53) and (2.70), 
we chose a simple, two- segment, acoustic and oceanographic model. The model was 
generated by dividing the acoustic path from the Northeast mooring to the Southeast 
mooring into two, flat-bottomed segments. Segment number 1 is north of the BSPF, 
while segment number 2 is south of the front. Each segment’s depth is the average 
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Figure 5-1: Bathymetry of two segment range dependent model from NE array to SE 
array. Solid line is the depth of the model. Dashed line is the depth at each CTD 
station. Segment #1 goes from 0-26. 93km. 

of the depths recorded at each of the CTD stations between the moorings. Figure 

5.1 shows the bathymetry (using CTD station depths) superimposed on the segment 
depths. The physical properties (temperature, salinity, et al.) are the averages from 
CTD stations between the moorings. Figures 5-2(a) and (b) show the average sound 
profiles of each segment. 

The total travel time disturbance will be taken as the sum of the disturbances in 
the two segments. This is somewhat inconsistent in that it ignores mode coupling 
effects due to the front. However, the exact magnitude of the mode coupling due to 
the front is not well known at this time, so we did not pursue this correction to our 
result. 

5.1.1 Acoustic Ray Results 

The simplest, shallow ocean model, given a strong mixed layer, is a two-layer ocean. 
In this model, the sound speed of each layer is the vertical average of the Fig. 5-2 
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Figure 5-2: (a) Averaged sound speed profile for segment # 1, and (b) Averaged 
sound speed profile for segment # 2 



sound speed profile for a given range segment. Table 5.1 lists the details for each 
segment. 

The travel time perturbations in ray theory can easily be calculated by perturbing 
the layer an amount equal to the peak to peak particle displacement £ and inputting 
the resultant sound speed perturbation SC into a discretized form of equation (2.7). 
This gives 



8T = £ 2SC(z n )NAz 



(5.1) 



sin x(z n ) 

where N is the number of complete ray cycles, and \ is the grazing angle. 

This calculation was run on segment # 1 for the M2 internal tide using a peak- 
to-peak particle displacement £ of 11 m with a source and receiver in the lower layer. 
It was assumed to first order that neither \ nor N changed, that is the ray path was 
identical before and after the sound speed profile was perturbed. Figure 5-3(a) shows 
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the results for eigenrays which penetrate the upper (mixed) layer. 

Figure 5-3(b) is the calculation performed for segment # 2 with a peak-to-peak 
particle displacement £ of 15 m, which was the maximum observed displacement of 
the thermocline seen in Fig. 4-3. The displacement of the thermocline due to the 
M 2 internal tide may be as small as 4 m, which was the computed M2 internal tide 
amplitude at the SW array using Fourier analysis. Using the smaller value of 4 m 
would scale Fig. 5-3(b) by one fourth. 

The calculation can also be performed using a more sophisticated ray tracing 
program such as MPP with fully range dependent profiles and geometry. As we are 
only interested in order of magnitude estimates with this ray calculation, we have not 
pursued this. 

5.1.2 Acoustic Normal Mode Results 

In calculating the acoustic modal results, the same two-segment view of the acoustic 
path was taken. Figures 5-4(a) and (b) show buoyancy frequency profiles calculated 
from equation (2.25) using averaged temperature and salinity profiles for segments # 
1 and #2. 

Using these buoyancy frequency profiles, equation (2.27) can be solved to give the 
internal tide and internal wave modes. Figures 5-5(a) and (b) show the first three 
internal tide modes for each segment at the M2 frequency. Similarly, given the sound 
speed profiles, equation (2.41) gives the adiabatic normal modes (Figs. 5-6(a) and 
(b)). In comparing Fig. 5-5 with Fig. 5-6, the case can be made that the third M2 
internal tide mode will possibly have the greatest effect on the acoustics, because its 



Table 5.1: Two-layer, range independent model 
Segment Nr. 1 Segment Nr. 2 



Distance (m) 


26,930 


7680 


Avg. water depth (m) 


185 


260 


Mixed layer depth (m) 


25 


25 


Sound speed (m/s) upper layer 


1473 


1479.9 


Sound speed (m/s) lower layer 


1453 


1465.1 
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Travel time perturbation (msec) Travel time perturbation (msec) 





Figure 5-3: Travel time perturbation of eigenrays versus launch angle. 
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Figure 5-4: Smoothed buoyancy frequency N profiles for (a) segment # 1 and (b) 
segment # 2 





- 0.1 - 0.05 0 0.05 0.1 

Amplitude 



- 0.1 - 0.05 0 0.05 0.1 

Amplitude 



Figure 5-5: First three internal tide modes for (a) segment # 1 and (b) segment # 2 
at uj = M2 
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Amplitude 



Amplitude 



Figure 5-6: First three adiabatic normal modes for (a) segment # 1 and (b) segment 



# 2 . 



maximum occurs near the maximum of the acoustic modes. 

The results of calculations using equation (2.53) are shown in Figs. 5-7 to 5-9 
in three different ways. First, as shown in Fig. 5-7, we assume that the lowest ten 
internal tide modes have an equal amplitude of 5.2 m. This calculation allows us 
to see how the internal tide and acoustic modes will interact, everything else being 
equal. For segment #1, the maximum travel time perturbation, indeed, occurs for 
the third internal tide mode. Except for the first internal tide mode, low internal 
tide modes dominate the scattering. The internal tide modes have little effect on 
the first acoustic mode. For segment #2, the low internal tide modes, including the 
first internal tide mode, also dominate the scattering, but the first acoustic mode is 
somewhat more affected by the internal tides. It should be noted that one reason 
that the travel time perturbations for segment # 2 are much smaller than those for 
segment # 1 is that the range for the former is four times smaller than the latter. 

For the second method, shown in Fig. 5-8, only the three lowest internal tide 
modes are used, instead of using the lowest ten modes as before. The three internal 
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Figure 5-7: Travel time perturbations at M2 using first ten internal tide modes 
equally weighted with a maximum amplitude of 5.2m. Lines are acoustic modes: 
#1(°), #2(+), #3(x), #4(»), #5(.). 
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tide modes are no longer equally weighted, but are now weighted by the amplitude 
factors, 0.58, 0.26, and 0.16, respectively, which were calculated in Section 4.2. This 
calculation allows us to see the contribution of each internal tide mode to the acoustic 
travel time perturbation. In segment #1, we see again that first internal tide mode 
has little effect on the low acoustic modes, and that the low internal tide modes (mode 
two and three) dominate the low acoustic modes. In segment #2 the first internal 
tide mode is dominant for all the acoustic modes calculated. 

For the third method, the three lowest internal tide modes are coherently summed 
together using the signs calculated in Section 4.2. The result is the solid line in Fig. 
5-9. The dashed line in Fig. 5-9 is the contribution from the first internal tide mode 
(the solid line in Fig. 5-8). In both segments #1 and #2, the first mode is the 
dominant contributor to the travel time perturbation. For segment #1 (Fig. 5-8) the 
magnitudes of the second and third internal tide modes are greater than that of the 
first internal tide mode over the range of the first four acoustic modes, but because 
the signs of the second and third internal tide modes are opposite at depths below 60 
m (where the group velocity dispersion is a maximum at the lower acoustic modes), 
the effects of the second and third internal tide modes largely cancel each other out, 
thus making the first internal tide mode the largest component of the travel time 
perturbation. For segment #2 this effect is less pronounced, because the magnitudes 
of the second and third internal tide modes are dissimilar. Figure 5-10 gives the 
root mean square of the travel time perturbations, the square root of the solution to 
equation (2.70), for segments #1 and #2 using the first three internal tide modes. 
The rms of the variance of the travel time perturbations for both segments is only 
slightly smaller than the respective travel time perturbations. 

5.1.3 Comparison of Ray and Modal Results 

For segment # 1 the Fig. 5-9(a) modal travel time perturbations of 0.5 to 35.5 msecs 
are of the same order as the ray travel time perturbations of 7.6 to 21.3 msecs. For 
segment # 2 the Fig. 5-9(b) modal travel time perturbations of 0.7 to 6.4 msecs are of 
the same order as the ray travel time perturbations of 2.7 to 4.1 msecs. This gives us 
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Travel time perturbation (msec) Travel time perturbation (msec) 





Figure 5-8: Travel time perturbations at M2. The first 3 internal tide modes are 
shown and are weighted by 0.58, 0.26, and 0.16, respectively. 
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Figure 5-9: Travel time perturbations at M2 using the sum of the first 3 internal tide 
modes, weighted by 0.58, 0.26, and 0.16. Shown for comparison (as a dashed line) is 
the contribution from the first internal tide mode at M 2 (Fig. 5-8). 
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Figure 5-10: RMS of the travel time perturbations at M2 using the first three internal 
tide modes. 
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some confidence that our calculations are at least consistent, within the assumptions 
we have made. 



5.2 Acoustic Normal Mode Results for the Inter- 
nal Tide and Internal Wave Fields 

Section 5.1.2 evaluated the travel time perturbations and covariances due to the 
largest internal tide component, the A/2 tide. In this section, using the same two 
segment model we perform the same calculations over the internal tide and internal 
wave band. Figures 5-1 1(a) and (b) are the travel time perturbations over the entire 
IW band for segments # 1 and #2. Although the travel time perturbation over the 
entire internal wave band has a similar shape to the travel time perturbation from 
the first IW mode at the A/2 tidal frequency (dashed line in Fig. 5-11), the internal 
wave band contributes significantly to the travel time perturbation. Modal travel 
time covariances were then computed using equation (2.70). The lowest ten acoustic 
modes, and the lowest 3 internal tide and internal wave modes were used in the 
calculations. The equation was discretized and summed from the inertial frequency 
/ to the buoyancy frequency N. The calculation is a conservative estimate, in that 
in the internal tide band, the tidal components were treated as delta functions, and 
only the A/2 tidal component was included. The 4 cpd and 8 cpd components were 
omitted, as their significance is doubtful (see Fig. 4-9). The internal wave band 
from 10 cpd to the buoyancy frequency N was treated as a continuous spectrum, and 
numerically integrated. Figure 5-12 gives the results for the acoustic modes. The rms 
of the variance of the travel time perturbations for both segments has only increased 
slightly, when the higher frequencies in the internal wave band are included. This is 
because the internal wave band is an order of magnitude down compared to the A/2 
tide in spectral amplitude. 

When we look at directional effects on the internal waves and tides, the dominance 
of the internal tide signal over the internal wave signal becomes even greater. At the 
A/2 tidal frequency, scattering due to waves in the across shelf direction ( 6 = 0) is only 
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Figure 5-11: Travel time perturbations using the first three IW modes from / to N. 
They are weighted by 0.58, 0.26, and 0.16. Shown for comparison (as a dashed line) 
is the contribution from the first internal tide mode at M2 (Fig. 5-8). 



65 





Figure 5-12: RMS of the travel time perturbations, using the first three IW modes, 
weighted by 0.58, 0.26, and 0.16, and from / to N. The dashed line is the rms of 
the travel time perturbation for the M2 component. The dotted line is rms of the 
travel time perturbation for the across shelf case. For segment #2 the dotted line lies 
directly over the dashed line. 
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slightly less than by waves going along shelf (0 = 90), since for the long wavelength 
of the first M2 mode (260km), the sine 2 function has a value of 0.967 as opposed to 
1.0 along shelf. (Thus internal wave scattering is roughly independent of 0\) For the 
internal waves, looking at 10 cpd as a representative frequency, the across shelf sine 2 
function is 0.0001 (compared to 1.0 for along shelf) due to the higher frequency and 
smaller wavelength. This compounds the order of magnitude difference in spectral 
amplitudes. 



5.3 Movement of the Front 

It was pointed out earlier in Section 3.1 that the BSPF can move up to 10km over a 
tidal cycle. To first order, this change can be approximated by taking the differences 
of the average sound speed between each segment, acting over one kilometer. The 
result is a travel time perturbation of 4.9 msecs/km, or up to 49 msecs over a tidal 
cycle. As the acoustic array is not expected to be able to filter out any modes higher 
than mode #5, the modal travel time perturbations calculated in Sections 5.1 and 
5.2 are not expected to be a significant source of noise for modal inverses for frontal 
structure. Internal wave noise could be significant, however, for the ray (higher mode) 
arrivals, a subject needing further investigation. Also, the phase relationship between 
the perturbation in travel time caused by the internal tide and internal wave fields 
and the movement of the front is unknown, because of the complex relationships 
involved in the bathymetry. Finally, due to the fact that one does have a fair degree 
of horizontal resolution in the cross front acoustic maps made tomographically, we 
may not be faced with a significant ambiguity between the vertical movement of 
the internal waves and the horizontal movement of the front. That is, the internal 
tide/internal wave “noise” may not significantly degrade the maps being made of the 
front. 



67 



Chapter 6 



Conclusions and Future Work 



As stated in the beginning, our goal in this thesis was to examine using acoustic to- 
mographic techniques to study internal waves and internal tides. Towards that end, 
we formulated an expression for the travel time perturbation of an acoustic wave by 
the internal tide and internal wave fields. The perturbation is maximum when the 
acoustic wave fields and internal wave fields are propagating at right angles to each 
other. Using temperature records at three deployed moorings, along with the tem- 
perature profiles from CTD casts, particle displacement records were obtained. We 
then determined the internal tide and internal wave spectra using Fourier analysis 
techniques. The shape of the particle displacement amplitudes at the SW mooring 
allowed us to estimate the amplitude weights of the first three internal wave modes. 
Knowing the amplitude spectra we then predicted the acoustic travel time pertur- 
bations. We found that for internal wave modes of equal amplitude that the third 
internal wave mode has the most effect on the acoustic modes. For the M2 tide, the 
first internal wave mode caused most of the travel time perturbation and variance. 
When the internal wave band was added, little change from the M2 result was seen 
for the first five acoustic modes. Since the acoustic VLA is not expected to be able 
to resolve modes higher than the fifth acoustic mode, for this application the use of 
the first internal wave mode at the M2 tidal frequency contains the bulk of the travel 
time perturbation and variance. We next compared the magnitudes of the travel time 
perturbation to the estimated change in the front due its movement. We concluded 
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that since the front could move upwards of 10 km during a tidal cycle, which would 
result in a travel time perturbation of up to 49 msecs, that the travel time perturba- 
tions due to the internal tide and internal wave fields would be a small, second order 
effect. 

Whereas in this thesis we have used Fourier analysis techniques to extract the 
internal tide and internal wave amplitudes, there exist several other methods which 
could be employed with varying degrees of success. One method is the Lagrangian 
contour following method used by Wang et al. [1991] that traces the minimum tem- 
perature through vertically interpolated temperature records. Another method is to 
determine the vertical velocities using the principle of the conservation of temperature 
which Holloway [1984] gives as 



dT/dt 

dTjdz 



(6.1) 



Both methods are more suited to situations where there are several temperature 
records in the vertical and where a minimal amount of smoothing needs to be done 
to remove the “spiking” due to eddies and other oceanographic processes. There also 
exist more sophisticated methods of tidal harmonic analysis which can be explored 
to extract the amplitudes such as those used by Foreman [1977] (cited in Holloway , 
1984). Finally, a last method to estimate the height of the M2 internal tide would 
be to calculate the movement of the thermocline using the anchor station CTD casts 
at the SW array. 

The analysis was hindered by our ignorance of the direction of propagation of the 
internal wave field. Levine and Richman [1989] have devised a method, whereby using 
a single mooring of current meters and thermistors, the direction could be resolved. 
Perhaps such an analysis of both the current meter and the temperature information 
may result in resolving the propagation direction. This information would allow us 
to check our calculations versus the measured data with much less uncertainty. 

Much more work remains to be done in fully understanding the nature of the 
internal tide and internal wave fields. The internal wave spectrum could be compared 
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to the Garrett-Munk internal wave model. The travel time perturbation covariance 
was computed for the 0 = 0 and 9 = 90 cases, but could also be computed for the 
isotropic case. Finally, future work in the Barents Sea could be done to resolve the 
phase relationship of the the displacement of the BSPF to the barotropic tide and 
to the internal tide, thereby allowing a more accurate determination of the acoustic 
signal travel time perturbation. 
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